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A correction for the Hartree-Fock Density of States for Jellium without Screening 
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We revisit the Hartree-Fock (HF) calculation for the uniform electron gas, or jellium model, whose 
predictions - divergent derivative of the energy dispersion relation and vanishing density of states 
(DOS) at the Fermi level - are in qualitative disagreement with experimental evidence for simple 
metals. Currently, this qualitative failure is attributed to the lack of screening in the HF equations. 

Employing Slater’s hyper-Hartree-Fock (HHF) equations, derived variationally, to study the 
ground state and the excited states of jellium, we find that the divergent derivative of the energy 
dispersion relation and the zero in the DOS are still present, but shifted from the Fermi wavevector 
and energy of jellium to the boundary between the set of variationally optimised and unoptimised 
HHF orbitals. The location of this boundary is not fixed, but it can be chosen to lie at arbitrarily 
high values of wavevector and energy, well clear from the Fermi level of jellium. 

We conclude that, rather than the lack of screening in the HF equations, the well-known qualitative 
failure of the ground-state HF approximation is an artifact of its nonlocal exchange operator. Other 
similar artifacts of the HF nonlocal exchange operator, not associated with the lack of electronic 
correlation, are known in the literature. 


INTRODUCTION 


The uniform electron gas, or jellium model, is an 
archetypal example in solid-state physics and many-body 
theory. Its treatment, in the Hartree Fock (HF) approxi¬ 
mation, can be found in classic textbooks [1-6], where, we 
learn that the HF equations applied to the ground state 
of the jellium, admit plane wave solutions with energy- 
wavevector dispersion relation given by, 
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kp is the Fermi wavevector, kp = 3n 2 (N/V). The single¬ 
particle energy e(k) is the sum of the free-electron energy, 
k 2 / 2, and the single-particle exchange energy. The Fermi 
wavevector kp is often expressed in terms of the mean 
radius per particle r s = y/9ir/Akp [5]; for typical values 
of r s in metals, the two terms in (1) are comparable in 
size. 

It is well known in the literature that the dispersion 
relation (1) has a logarithmically divergent derivative at 
the Fermi energy, shown in Fig. 1. Another marked dif¬ 
ference between the free electron result and the HF so¬ 
lution for jellium, evident in Fig. 1, is the considerably 
increased bandwidth of the HF dispersion. Finally, it 
is well known that in the HF approximation the DOS 
for jellium vanishes at the Fermi level (Fig. 1), since the 
DOS is inversely proportional to the derivative of the dis¬ 
persion. The zero in the DOS at the Fermi level suggests 
that jellium is a semimetal, in obvious disagreement with 
experimental evidence for simple metals, such at sodium 
or aluminium, which are described accurately by the jel¬ 
lium model. 

In the literature, the qualitatively wrong description of 
jellium in the HF approximation is attributed to the long 
range of the Coulomb repulsion [1-6]. It is well known 
that the flawed description can be corrected by the intro- 




Figure 1: Solid lines show ground-state HF results for 
jellium, compared to free-electron results in dotted 
lines. (r s /ao = 4, = kp/ 2.) Top: Energy vs wave 

vector dispersion relation e(k). The logarithmic 
divergence in the derivative, de/dk , is marked with a 
triangle (A). Bottom: DOS, showing the unphysical 
zero at the Fermi level for jellium. 


duction of electronic many body correlation effects [1 7], 
which screen the bare Coulomb potential and thus elimi¬ 
nate the unphysical divergent derivative of the dispersion 
relation, the zero in the DOS at the Fermi level, and also 
reduce the bandwidth of the HF dispersion relation of 
jellium. 

In an effort to understand whether HF’s lack of screen¬ 
ing actually plays a role, we revisit the HF study of jel- 
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lium, attempting to correct the qualitative errors of the 
HF description, but without including any form of elec¬ 
tronic correlation. For this purpose, we employ Slater’s 
hyper-HF (HHF) theory for the ground and the excited 
states of an IV-electron system [8]. Specifically, we use 
the single-particle HHF equations by Gidopoulos and 
Theophilou [9, 10], who considered an IV-electron system 
described by a Hamiltonian H and then variationally op¬ 
timised the average energy ]T n (<i> ra |H|<I> ra ) °f configu¬ 
rations (IV-electron Slater determinants $„) constructed 
from a basis set of R spin-orbitals, R> N. 


THE HYPER-HARTREE-FOCK EQUATIONS 
FOR JELLIUM 

The aim in HHF theory is to obtain approximations, 
at the HF level of description, for the ground and the 
excited states of an IV-electron system. These states 
are represented by IV-electron Slater determinants, con¬ 
structed from a common set of spin-orbitals. Obvi¬ 
ously, to have the flexibility to describe excited states, 
the number of spin-orbitals ( R ) must exceed the num¬ 
ber of electrons. For example, say we are interested to 
approximate the ground and excited states of the he¬ 
lium atom. In the HF ground state of He, the Is orbital 
{(pis) is doubly occupied. To study a couple of excited 
states, we need at least one more spin-orbital and the 
next one is: ip\ s . With the three available spin-orbitals, 
IsT, Is-*-, 2^, we can form three configurations for the He 
atom (two-electron Slater determinants): = [Is - * - , Is 1 "], 

$2 = [ls^,2st], $ 3 = [lst,2sl]. In HHF theory, we vari¬ 
ationally optimise the three common spin-orbitals simul¬ 
taneously, by minimising the sum of the expectation val¬ 
ues The minimisation leads to the HHF 

single-particle equations for the three spin-orbitals. It 
turns out that these equations resemble the ground-state 
HF equations for the lithium atom (three electrons) but 
with a weakened Coulomb repulsion between the three 
electrons, to keep balance with the nuclear charge which 
has remained that of the He nucleus. 

In general, in HHF theory [8, 9] for an IV-electron sys¬ 
tem, one considers a set of R orthonormal spin-orbitals, 
with R > N. On this spin-orbital basis set one may 
define, D = R\/(N\(R — IV)!), IV-electron Slater deter¬ 
minants. 

The derivation of the single-particle HHF equations in 
Ref. 9 is based on Theophilou’s variational principle [11], 
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where {E%} are the D lowest eigenvalues of the IV- 
electron Hamiltonian H. 

An extension of the variational principle, with unequal 
weights in the sums in (2) was proposed by Theophilou 


[12], and independently by Gross, Oliveira, Kohn [13]. 
These variational principles can be derived as special 
cases from the Helmholtz variational principle in statis¬ 
tical mechanics [14, 15]. In particular, the inequality in 
(2) arises as the high temperature limit of the Helmholtz 
variational principle. 

Optimisation of the R spin-orbitals {yft} to minimise 
the sum of the energies on the l.h.s. of (2) leads to the 
following single-particle equations for the spatial part of 
the spin-orbitals [9] (in atomic units): 
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are the Coulomb and exchange operators respectively. 
Kxt signifies the attractive potential of the nuclear 
charge. For R = IV, Eqs. 3 reduce to the familiar ground- 
state HF equations. In Eqs. 3, the orbitals ipi, with 
i = 1,...,J?, are correctly repelled electrostatically by 
a charge of IV — 1 electrons. In contrast to the ground- 
state HF case [16], this holds true even for the orbitals 
that are not occupied in the HHF ground-state Slater de¬ 
terminant as long as these orbitals are variationally op¬ 
timised, i.e., for ipi, with N < i < R. Furthermore, the 
orbitals that are left variationally unoptimised, i.e., ipi , 
with i > R, are repelled by a charge of N — 1 + (1/A) elec¬ 
trons. In the HHF equations, the well-known asymmetry 
in the treatment of the variationally optimised and unop¬ 
timised orbitals by the nonlocal exchange operator [16] 
is still present, but softened (for large A), compared with 
ground-state HF. We note that for R > IV, Koopmans’ 
theorem [19, 20] ceases to hold for the HHF equations. 

The HHF equations (3) have the form of ground-state 
HF equations for a virtual system of R-electrons, where 
the electronic Coulomb repulsion is multiplied by 1/A: 
1/|r — r'| —> A _1 /|r — r'|. Therefore, the calculation 
of the optimal spin-orbitals to represent the ground and 
excited-states of an IV-electron system in the HHF ap¬ 
proximation, is reduced to the calculation of the ground- 
state HF orbitals of a fictitious system with a greater 
number of electrons R > N, and scaled down electronic 
Coulomb repulsion. A related approach is the “super- 
hamiltonian method” by Katriel [22, 23]. 
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Finally, before applying the HHF equations to jellium, 
we remark that correlated, approximate eigenstates of 
the Hamiltonian H can be obtained by diagonalising 
the matrix |$ m ) [9], where $„ are the ./V-electron 

HHF Slater determinants. This configuration-interaction 
method employs the HHF spin-orbitals, which are opti¬ 
mised to represent on equal footing the ground and the 
excited states of H. 


Solution of the HHF equations for jellium. 

Similarly to the HF ground state, the HHF equations 
for jellium admit plane wave solutions. It follows that the 
ground-state ./V-electron Slater determinant and its total 
energy are the same in the HF and HHF approximations. 

Although the HF and HHF equations for jellium admit 
the same solution for the orbitals, the dispersion relations 
for the single-particle energies e{k) and A (k) differ. In 
particular, the HHF dispersion, A (k), results from an op¬ 
timisation that involves a broader range of wavevectors 
than the HF dispersion e(k). 

Following the standard treatment in textbooks [1-6], 
it is straightforward to work out directly the solution of 
the HHF Eqs. 3. Here, we exploit the similarity of Eqs. 3 
with ground-state HF equations of an A-electron system, 
to obtain that the HHF dispersion relation, A (k), will 
be given by (1) with the single-particle exchange energy 
scaled down by the factor 1/A: 
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kR is the Fermi wave vector of the virtual I?-electron 
system, 

4 = 3tt 2 ^ . ( 8 ) 

Dividing k F /k F . and taking the thermodynamic limit, 
N, V —> oo, with the ratio A fixed, we obtain: 


kR = A l ^k F ■ 


(9) 


Substitution of the above into Eq. (7) yields the desired 
expression for the single particle energy levels of jellium, 
in terms of A and the Fermi wavevector kF of the actual 
./V-electron system: 
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The dispersion relation in Eq. (10) reduces to the 
ground-state Hartree-Fock result for A = 1, and to the 
free electron dispersion relation, A(fc) = k 2 / 2, in the limit 


A —> oo (See Fig. 2). For increasing A, the bandwidth 
of the HHF dispersion, A (k), decreases compared to the 
ground-state HF dispersion, e(k). For A —>• oo the ex¬ 
change term in HHF dispersion vanishes and A (k) reduces 
to the free-electron result. 

Importantly, the wave vector at which the logarithmi¬ 
cally divergent derivative occurs is shifted from k F to kR , 
such that the divergence no longer occurs at the Fermi 
energy of the physical ./V-electron system, when the num¬ 
ber of optimised orbitals is R > N. 
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Figure 2: Excited-state HF results for jellium, for 
various A = R/N, compared to free electron results 
(dotted lines). (r s /ao = 4, e° F = kip/2.) Top: Energy vs 
wave vector dispersion relations A (k). When A = 1, 
A(fc) = e(k). Triangles (A) mark logarithmic divergence 
in dX/dk at Fermi level of fictitious Abelectron system. 
Bottom: DOS g(X(k)) : showing the zero at the Fermi 
level of the fictitious I?-electron system. 


The DOS, g(X)SX, can be obtained directly from Eq. 10 
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[2] and is given by the parametric equation: 
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g(X(k)) is expressed in terms of k R (rather than A 1 / 3 ^) 
to keep the notation simple. For any finite A, the DOS 
still vanishes. However, as shown in Fig. 2, the zero in 
the DOS occurs at the Fermi energy of the fictitious R- 
electron system, A (k R ), rather than the Fernri energy of 
the physical system A (kp). 

DISCUSSION 

In metals, screening is an important effect that reduces 
the range of the effective repulsion between electrons, 
shielding any charge at distances greater than a charac¬ 
teristic screening length. In the literature of many-body 
theory [1, 3-6] and solid-state physics [2-4], where jellium 
is a paradigm, the qualitatively flawed description of met¬ 
als by the HF approximation is attributed to the long- 
range nature of the Coulomb interaction, which, com¬ 
bined with the neglect of correlation, deprives from the 
HF equations the flexibility to model the phenomenon of 
screening. This understanding of HF’s failure is further 
supported by the softening of the divergence in the slope 
of e(k), after replacing the bare Coulomb potential in 
the HF nonlocal exchange term by a screened Coulomb 
potential [2]. 

On the other hand, in the theoretical chemistry liter¬ 
ature, it is well known that the HF nonlocal exchange 
term, in finite systems, gives rise to several counterintu¬ 
itive results, reminiscent of the HF anomalies in jellium, 
which are not associated with HF’s lack of correlation. 
For example, Handy et. al [17] disproved the widely 
held view [1], that the HF nonlocal exchange potential 
decays as (— 1/r) at large distances. In particular, in 
Ref. 17 it was demonstrated that the asymptotic decay 
of an occupied HF orbital ifi with eigenvalue £j is not 
~ exp(— \J—2eir), as would be expected from the (—1/r) 
tail of the exchange potential. On the contrary, in HF 
all the occupied orbitals <fi decay uniformly at large dis¬ 
tances away from the system, regardless of their energy 
eigenvalue. 

As already mentioned, it is also widely known that 
the HF exchange operator deals with the occupied and 
unoccupied orbitals in the ground-state HF Slater de¬ 
terminant in an asymmetric way [16]: for an IV-electron 
system, the (variationally optimised) occupied orbitals 
are repelled electrostatically by a charge of N — 1 elec¬ 
trons, while the (variationally unoptimised) unoccupied 
orbitals feel the stronger repulsion of N electrons, mak¬ 


ing the unoccupied orbitals too diffuse, and raising their 
energy eigenvalue to unphysically high values [16]. 

It follows that the energy to excite an electron from an 
occupied HF orbital ifi to a virtual HF orbital tp a , keep¬ 
ing the other occupied orbitals frozen, must be smaller 
than the eigenvalue difference e a — £*. This is because 
the energy e a incorporates the Coulomb interaction of 
the orbital ip a , hosting the electron after excitation, with 
all the occupied orbitals, including the orbital ipi accom¬ 
modating the same electron before excitation. In this 
sense, the Coulomb interaction of the orbitals tpi and (p a 
can be interpreted as a form of self-interaction raising 
the energy of the virtual level s a . It is similar to the 
“ghost” self-interaction discussed in Ref. 21. This inter¬ 
pretation is consistent with viewing the virtual HF en¬ 
ergies as single-particle levels of the iV-electron system. 
Alternatively, by extending the proof of Koopmans’ the¬ 
orem [19], it can be shown that e Q is equal to the negative 
of the electron affinity of the system to bind an electron 
at the virtual orbital ip a , see e.g. Ref. 20. The interpreta¬ 
tion of the virtual energies as negative affinities amounts 
to regarding the unoccupied levels as virtual levels of an 
(.N + l)~electron system. 

We note that for jellium, where the occupied and vir¬ 
tual orbitals are plane waves, the self-interaction error 
of the virtual energy levels (discussed above) vanishes 
in the thermodynamic limit. Therefore, it makes sense 
to consider that the HF unoccupied levels represent vir¬ 
tual single-particle levels of the IV-electron system, and 
to study their dispersion relation and density of states. 

A consequence of the asymmetry in the treatment of 
occupied versus unoccupied orbitals in the HF equations, 
is presented by Bach et al. [18], who prove that, for 
a finite system, the highest occupied spin-orbital in the 
ground-state unrestricted-HF Slater determinant, is non¬ 
degenerate: a nonzero gap separates it from the lowest 
unoccupied spin-orbital, even in systems with an odd 
number of electrons, in contrast to physical expectation. 

Recently, Hollins et al., using two different meth¬ 
ods, the optimised effective potential, or exact exchange 
method [24], and the local Fock exchange potential 
(LFX) method [25], showed that it is possible to omit cor¬ 
relation and still obtain an accurate dispersion relation 
e{k) for simple metals, provided the exchange potential 
term is local v x = v x (r) [25]. The LFX potential, defined 
as the local exchange potential with the HF ground-state 
density [25, 26], is particularly interesting in our context: 
even though the determination of the LFX potential is 
based on the same ground-state HF calculation that gives 
a very poor prediction for the dispersion of simple metals 
(e.g. Na, Al), the band-structure of the LFX potential 
[25] almost coincides with the band-structure of the local 
density approximation in density functional theory [27], 
which, by construction, is very accurate for these sys¬ 
tems. With regard to the vanishing of the HF DOS of 
metals at the Fermi energy, Hollins et al. [25] argue that 
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it appears to be an artifact of the HF nonlocal exchange 
operator. 

Our work supports this point of view, by studying the 
nonlocal HF exchange term directly. We find that the 
well-known anomalies of the HF description of jellium are 
still present in the solution of the HHF equations. How¬ 
ever, the location of the divergent derivative of A(fc) and 
the zero in the DOS are no longer at the Fermi level of the 
actual iV-electron system. Instead, they are positioned 
at the border separating the variationally optimised and 
unoptimised orbitals. By choosing to variationally opti¬ 
mise a very large number of orbitals, the unphysical zero 
in the DOS can be pushed to very high energies, avoid¬ 
ing completely the window of single-partice energies that 
can be of any relevance to the ground state of jellium. 
Therefore, it is no longer justified to relate the mobile di¬ 
vergence in the derivative of the HHF dispersion and the 
travelling zero in the HHF DOS of jellium, with the lack 
of electronic correlation in the HHF approximation, even 
though these anomalies can still be removed by screening 
the Coulomb repulsion in the exchange term. 

This strengthens the view that the failure of the 
ground-state HF solution for jellium is also an artifact 
of the nonlocal exchange operator. At least part of the 
explanation of the divergent derivative of the disper¬ 
sion seems to be the asymmetry in the treatment of the 
variationally optimised and unoptimised orbitals: as the 
wavevector k crosses kn from below, the plane wave solu¬ 
tions of the excited-state HHF equations are subjected to 
a discontinuous drop of the nonlocal exchange term and 
a correspondingly discontinuous increase in the Coulomb 
repulsion from the Hartree term. This raises the single¬ 
particle eigenvalue A (k) to higher energies, diminishing 
the DOS in the neighborhood of A(/cr). The same mech¬ 
anism operates in the divergent slope of the ground-state 
HF dispersion e{k). 

In conclusion, we suggest that the qualitative failure 
of the HF approximation for jellium is unrelated to the 
lack of correlation in the HF approximation. Instead, 
this failure is another example in the list of counterintu¬ 
itive results caused by the nonlocality of the HF exchange 
operator. Our work complements the work of Hollins et 
al. in Ref. 25 and our conclusion is contrary to what is 
currently written in almost any textbook in the fields of 
many body theory and solid state physics. To the best 
of our knowledge, Slater gave the only hint so far in the 
literature that the HF failure may be unrelated to elec¬ 
tronic correlation. In his textbook on “Insulators, Semi¬ 
conductors and Metals” [28] he writes that “it cannot be 
in accordance with experiment to write the total energy 
of the electronically excited state of the crystal as a sum 
of one-electron energies of the type of ...” (cf. the HF so¬ 
lutions). “A great deal of effort has gone into explaining 
this apparent paradox connected with the free-electron 
theory of electrons in metals. ... It is the impression of 
the present author, however, that we do not need to look 


for any deep and profound explanation”. 
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